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We study a driven-dissipative system of atoms in the presence of laser excitation to a Rydberg 
state and spontaneous emission. The atoms interact via the blockade effect, whereby an atom in 
the Rydberg state shifts the Rydberg level of neighboring atoms. We use mean-field theory to study 
how the Rydberg population varies in space. As the laser frequency changes, there is a continuous 
transition between the uniform and antiferromagnetic phases. The nonequilibrium nature also leads 
to a novel oscillatory phase and bistability between the uniform and antiferromagnetic phases. 



PACS numbers: 



The behavior of matter far from equihbrium is a fasci- 
nating area of study. The presence of driving and dissipa- 
tion can lead to remarkable phenomena that are not pos- 
sible in equilibrium. This has motivated much research 
on nonequilibrium physics in classical systems, such as 
fluids, chemical reactions, and biological media [l|, i- An 
interesting question is: what novel phases appear when a 
quantum system is driven far from equilibrium? Recent 
cold-atom experiments have studied equilibrium quan- 
tum systems in great detail, but they are also a natural 
setting to study nonequilibrium quantum systems due to 
the tunability of driving and dissipation [J-Q. 

In this paper, we study a nonequilibrium many-body 
quantum system interacting via Rydberg blockade. A 
Rydberg atom is one whose electron is excited to a 
high energy level n. The van der Waals interaction be- 
tween two atoms in identical Rydberg levels scales as 
n^-^, and this leads to a blockade effect for large n: when 
one atom is excited to the Rydberg state, it prevents 
nearby atoms from being excited. This is the basis for 
quantum information processing schemes with Rydberg 
atoms [gl-fl^ and a variety of novel phenomena (l5l - l22| . 
In these schemes, spontaneous emission should be min- 
imized, since it destroys quantum information. On the 
other hand, spontaneous emission as a source of dissipa- 
tion may lead to interesting physics, and it can actually 
be tuned by using different Rydberg levels. 

We study a lattice of atoms continuously excited to 
the Rydberg state and spontaneously decaying back to 
the ground state. Consider the Rydberg population of 
each atom, i.e., the fraction of time it spends in the 
Rydberg state. What is the spatial distribution of the 
Rydberg population in steady state? Using mean-field 
theory, we show that as the laser frequency is varied, 
the system undergoes a continuous transition between a 
phase with spatially uniform population and a phase with 
higher population on every other atom. We call the lat- 
ter the antiferromagnetic phase, since a two-level atom is 
formally equivalent to a spin-1/2 particle (ground and ex- 
cited states correspond to down and up spins) ^2^. The 
nonequilibrium nature also leads to a novel oscillatory 
phase, in which the Rydberg population oscillates peri- 
odically in time, and bistability between the uniform and 



antiferromagnetic phases. Simulations of the full quan- 
tum model in ID (where mean-field theory is least accu- 
rate) show that there are short-range antiferromagnetic 
correlations but not long range order. Our work can be 
extended to more general dipolar gases and NMR. 

First, we describe the Rydberg interaction [l3|- Sup- 
pose two atoms are in the same Rydberg state nlj. There 
is a dipole-dipole matrix element between \nljnlj) and 
nearby energy levels, and this interaction shifts the en- 
ergy of \nljnlj) by an amount V. When the atoms are 
separated by a small distance R, the dipolar interaction 
dominates {V ~ —C^/R^), but for large distances, the 
van der Waals interaction dominates {V « —Cq/R^). For 
mathematical convenience, we use the van der Waals in- 
teraction and a \nsi/2nsi/2) state, so that the interaction 
is short range and isotropic. However, it is straightfor- 
ward to extend the analysis to long-range and anisotropic 
interactions. The value of Cg depends on n,l,j and is 
tabulated in Refs. [M^. 

Consider a lattice of atoms that is uniformly excited 
by a laser from the ground state to a Rydberg state. The 
atoms are assumed to be fixed in space. Since the van 
der Waals interaction decreases rapidly with distance, we 
assume nearest-neighbor interactions. Let \g)j and \e)j 
denote the ground and Rydberg states of atom j. The 
Hamiltonian in the interaction picture and rotating- wave 
approximation is (ft = 1) 

H = |e)(e|,®|e)(e|fe, (1) 

H, = -A|e)(e|, + ^(|e)(.g|,-f |g)(e|,). (2) 

The second term in Eq. ^ is the Rydberg interaction, 
and Hj is the Hamiltonian for a two-level atom interact- 
ing with a laser. A = a;^ — cJq is the detuning between the 
laser and transition frequencies. 17 is the Rabi frequency, 
which depends on the laser intensity. 

The lifetime of the Rydberg state is limited by sev- 
eral processes: spontaneous emission, blackbody radia- 
tion, and superradiance [13] • We account for spontaneous 
emission from the Rydberg level using the linewidth 7. 
When a Rydberg atom spontaneously decays, it usually 



goes directly into the ground state or first to a low-lying 
state [2^; the low-lying states are relatively short-lived, 
so we ignore them. We also ignore blackbody radiation 
and superradiance, both of which transfer atoms in a 
Rydberg level to nearby levels. Blackbody radiation can 
be minimized by working at cryogenic temperatures psl j , 
and it is not clear if superradiance is important when the 
interaction V is large [27, 29]. Future treatments could 
account for them by considering several Rydberg levels 
instead of just one. 

Thus, each atom has two possible states, and the sys- 
tem is equivalent to a dissipative spin model. Previous 
works have added dissipation to other spin models by 
coupling each spin to a heat bath; in those works, there 
is global thermal equilibrium, and the spins are described 
by an effective partition function (sol. l3l|. However, in 
quantum optics, dissipation from spontaneous emission 
leads to a nonequilibrium situation, since the coupling 
to the heat bath is weak and Markovian [l^l- The den- 
sity matrix for the atoms, p, is described by a master 
equation that is local in time: 



(a) 



(b) 



p = -i[H,p]+C[p] 
1 
2 



{\e){e\j,p} + \g){e\j p\e){g\j 



(3) 
•(4) 



The nonequilibrium nature is exhibited in the interplay 
between unitary and dissipative dynamics 0, IE| j and we 
are interested in the properties of the steady-state solu- 
tion of Eq. ©. 

Due to the complexity of the full quantum problem, 
we use mean-field theory. For equilibrium spin mod- 
els, mean- field theory is useful for determining the ex- 
istence of different phases ^3,- Its predictions are ac- 
curate in high dimensions but not in low dimensions. 
For the current nonequilibrium case, we use the ap- 
proach of Refs. [1, factorize the density matrix by 
site, p = pj, and work with the reduced density ma- 
trices, pj = Tv^jp. This accounts for on-site quantum 
fluctuations but not inter-site fluctuations: for atom j, 
the interaction, |e)(e|j (g) J2k replaced with the 

mean field, |e)(e|j J2k Pk,ee- In high dimensions, this is a 
good approximation, since fluctuations of the neighbors 
average out. 

Then the evolution of each pj is given by 



Wj = — 2rilm Qj — jiwj + 1) , 



V 

qj = i A - — 2^(wfc + l) 



(5) 

I3 - . (6) 



where we have defined the inversion Wi = 



Pj,ee Pj,gg 



and 



off-diagonal element qj = Pj,eg- The Rydberg population 
Pj.ee — {wj + l)/2 is the observable measured in the ex- 
periment by measuring the scattering rate of each atom. 
Wj = —1 and 1 mean that the atom is in the ground and 
Rydberg states, respectively. Equations ([5]) and dH) are 
the Optical Bloch Equations, except that the Rydberg 
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FIG. 1: (a) When one atom is excited to the Rydberg state, 
it sliifts the transition frequency of a neighboring atom by V . 
(b) The lattice is divided into two sublattices. 



interaction introduces nonlinearity: the detuning for an 
atom is renormalized by the excitation of its neighbors 
(Fig.Uti). 

Since the system is dissipative, it will end up at an 
attracting solution, which can be a fixed point, limit cy- 
cle, quasiperiodic orbit, or strange attractor [3j]. (We 
have not observed the latter two.) We want to know: for 
given parameter values, how many steady-state solutions 
are there and are they stable? A solution is stable or 
unstable if a perturbation to it decays or grows, respec- 
tively; the system will end up only in a stable solution. 

Equations ([5]) and ([6|) always have a steady-state solu- 
tion, in which the Rydberg population is uniform across 



the lattice {wj 



qj = q). For some parameter values. 



this uniform solution is stable, but for others, it is unsta- 
ble to perturbations of wavelength 2. In the latter case, 
the lattice divides into two alternating sublattices, and 
the atoms on one sublattice have higher Rydberg popula- 
tion than the other. Hence an antiferromagnetic pattern 
emerges from the uniform solution through a dynamical 
instability. To simplify the discussion here, we keep track 
of only the two sublattices instead of every site (Fig.[TlD). 
We stress that the antiferromagnetic transition is not an 
artifact of using a bipartite lattice, as shown explicitly in 
the supplemental material 34 1. 

To simplify the equations, we rescale time by 7 and also 
rescale the Rabi frequency ft — ^/'j, detuning A = A/7, 
and interaction c = dyf-f = —dCe/'jR^, where d is the 
lattice dimension. Labeling the sublattices 1 and 2, 



Wi 

qi 

92 



—2Q Im qi — wi — 1 
—2Q Im (72 — ^2 ~ 1 

i [A - c(u;2 + 1)] qi - 



qi 



(7) 

(8) 

(9) 



l[A-ciwi + l)]q2~^+l-W2. (10) 



There are six nonlinear differential equations (since qi 
and 52 are complex) and three parameters (i7. A, c). The 
uniform version of these equations (wi — W2, qi = 92) 
has been studied before in the context of a medium that 
interacts with its own electromagnetic field; it is known 
that there is bistability [s^ . We are considering the more 
general case by letting the sublattices differ. 

In the supplemental material [s^ . we determine the 
solutions and stabilities for Eqs. ©-([TU]). Here, we sum- 
marize the main results. Consider first the fixed points. 
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i.e., when wi = W2 = qi = q2 = 0. There are two types 
of fixed points: the uniform fixed points (wi = ^2) cor- 
respond to spatially homogeneous Rydberg excitation, 
while the nonuniform fixed points (wi 7^ W2) correspond 
to the antiferromagnetic phase, i.e., when one sublattice 
has higher excitation than the other. 

There are either one or three uniform fixed points, cor- 
responding to the real roots of a cubic polynomial. 



(a) 



f{w) = C^W^ 



c(2A - 3c)w'^ 



2 



1 

4 



+ (A - 3c)(A -c)\w + {A-cf + -. (11) 

As the parameters change, pairs of uniform fixed points 
appear and disappear via saddle-node bifurcations. The 
uniform fixed points never undergo Hopf bifurcations, so 
we do not expect limit cycles emerging from them [33| . 

There are up to two nonuniform fixed points, given by 
the real roots of a quadratic polynomial. 



g{w) = 0^(1 + 4A2 + 217^)^2 

-2c[(A-c)(1+4A2)- 
+0^(1 +4A2) - 2cA(l 

-4A2 + 2172)2 _ 



(2A- c)rj2]w 
-4A2 + 21^2) 



4" 



(12) 



The two roots correspond to wi and W2- As the parame- 
ters change, the two nonuniform fixed points appear and 
disappear together. 

Since the laser detuning A is the easiest parameter 
to vary experimentally, we describe what happens as a 
function of it (Fig. [5]). Suppose A starts out large and 
negative. There is one stable uniform fixed point and 
no other fixed points. As A increases, the uniform fixed 
point may undergo a pitchfork bifurcation, in which it be- 
comes unstable and the nonuniform fixed points appear. 
The bifurcation is supercritical, which means that when 
the nonuniform fixed points appear, they are stable and 
coincide with the uniform fixed point [3j|. Thus, this is a 
continuous phase transition between the uniform and an- 
tiferromagnetic phases. As A increases further, there is 
another supercritical pitchfork bifurcation, in which the 
same uniform fixed point becomes stable again and the 
nonuniform fixed points disappear. As A increases fur- 
ther towards 00, there is again one stable uniform fixed 
point and no other fixed points. 

Although the nonuniform fixed points are stable when 
they appear and disappear, they could become unstable 
in between via a Hopf bifurcation ^33] . We find numeri- 
cally that sometimes the nonuniform fixed points do have 
Hopf bifurcations (Fig. ^^d) and give rise to a stable limit 
cycle, in which wi and •W2 oscillate periodically in time 
(Fig. [5^). This oscillatory phase is due to the nonequi- 
librium nature of the system. 

Thus, in mean-field theory, there are three phases: uni- 
form, antiferromagnetic, and oscillatory. FigurelJb shows 
a phase diagram in A, D, space. For some parameters. 





FIG. 2: Bifurcation diagram showing fixed-point solutions as 
function of A, with c = 5 and (a) fl = 0.5 and (b) Q = 1.5. 
The inversion w is -1 (1) when the atom is in the ground 
(Rydberg) state. Solid (dashed) lines denote stable (unstable) 
fixed points. Black (red) lines denote uniform (nonuniform) 
fixed points. Green points denote bifurcations. In (b), the 
nonuniform fixed points undergo Hopf bifurcations at A = 
3.48 and 1.33, and there is a stable limit cycle in that interval 
(shown in Fig. [3t-). 
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FIG. 3: (a) OsciUatory steady-state solution (limit cycle) for 
c = 5, fl = 1.5, and A = 1.5. (b) Phase diagram for mean- 
field theory in Q,A space for c — 5. The system is either in 
the uniform, antiferromagnetic, or oscillatory phase. It can 
be bistable between uniform and antiferromagnetic phases or 
between uniform and oscillatory phases. 



the system is bistable between uniform and antiferromag- 
netic or between uniform and oscillatory (Fig. [^b); the 
final state depends on the initial conditions. 

We also numerically solve the original master equa- 
tion, Eq. ([U, in ID, where mean-field theory is least ac- 
curate. We use fourth-order Runge-Kutta integration to 
find the steady-state p for a chain of length N — 10. Fig- 
ure 2^ shows the correlation as a function of distance, 
{EiE,+j) - {Ei){Ei+j), where Ei = |e)(e|i. The rapid 
decay suggests that there is no long range order in ID, 
but the fact that it alternates sign means that there is 
an antiferromagnetic tendency. We also calculate the 
order parameter, [{{Ee — Eo)^)]^^'^, where the operator 
= Gvon mcasurcs the average Rydberg pop- 

ulation on the even sublattice, and Eq does likewise for 
the odd sublattice. The order parameter measures the 
difference between the two sublattices: it is when they 
are identical (uniform phase) , but positive when they are 
different (antiferromagnetic and oscillatory phases) . The 
order parameter is largest for roughly the same parameter 
space, for which mean-field theory predicts the uniform 
phase to be unstable (compare Fig. with Fig. [SJd). 
Thus, mean-field theory captures some qualitative as- 
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FIG. 4: Numerical solution of master equation for ID chain 
of length A'^ = 10 with periodic boundary conditions. Steady 
state p is found after integrating for time 'yt = 20. Parameters 
are Q. — 1.5 and V — 67, which is equivalent to Fig. [Jla. (a) 
Correlation as a function of distance j for A = 0. (b) Order 
parameter as a function of detuning. 

pects of the full quantum model in ID, but it remains 
to be seen whether there is long-range order in higher 
dimensions, where mean-field theory is more accurate. 
Also, the prediction of an oscillatory phase suggests that 
the emitted light has system- wide temporal correlations; 
indeed, we have found strong photon correlations, which 
will be reported in detail elsewhere (36| . 

Since it is difficult to simulate large systems, exper- 
iments with atoms in an optical lattice could provide 
much information. For example, one can use ^""Rb and 
a two-photon excitation scheme to go from the ground 
state 5si/2 to the Rydberg state 23si/2, which has van 
der Waals interaction Cg = —870 kHz iim^ [2^ and 



linewidth j/2t: = 14.7 kHz at K A d-dimensional 
lattice with spacing i? = 1.5 /xm has interaction strength 
V = 76 kHz and c = 5. 2d. The Rydberg population 
of each atom may be measured by imaging the sponta- 
neously emitted photons; in the antiferromagnetic phase, 
every other atom fluoresces more. Alternatively, the 
ground-state population may be measured using repeated 
projective measurements on a 5s — 5p transition. A prac- 
tical setup would be to use a microscope that both pro- 
duces the lattice and images the atoms [37| . 

Thus, a driven-dissipative system of Rydberg atoms 
has a unique type of antiferromagnetism. The next step 
is to investigate in more detail how the full quantum 
model behaves in low dimensions. Our work can be ex- 
tended to Rydberg states with anisotropic and long-range 
interactions. Such interactions usually give rise to very 
rich physics [38l |. so the nonequilibrium version should 
be interesting. One can also see what happens when the 
atoms are not fixed on a lattice but free to move; this is 
reminiscent of classical reaction-diffusion systems [l], [1] . 
Finally, we note that a system of interacting Rydberg 
atoms is similar to a system of spins interacting with 
each other's magnetic dipolar field [H, H^. Thus, when 
an NMR system is made nonequilibrium with continuous 
driving and spin relaxation, the spins may form a stable 
pattern in space. 

We thank A. Daley and G. Refael for useful discussions. 
This work was supported by NSF Grant No. DMR- 
1003337. 
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SUPPLEMENTAL MATERIAL 

In this supplement, we provide details on the steady-state solutions of the mean-field model. In Appendix A, we 
study the solutions of the bipartite lattice. In Appendix B, we study the solutions of the complete lattice. We use 
extensively the language of nonlinear dynamics, which is explained in Ref. 



Appendix A: Bipartite lattice 

Here, the lattice is bipartite, and we only keep track of the two sublattices, which are labelled 1 and 2. The system 
is described by 

wi = -2r2Im (7i - wi - 1 , (Al) 

W2 = -2film - W2 - 1 , (A2) 

91 = ^ [A - c{w2 + 1)] 91 - Y + i^wi , (A3) 

q2 = i[A- c{wi + 1)] ^2 - Y + i^w2 , . (A4) 

Remember that wi , W2 are real while 91 , 92 are complex. There are six differential equations and three pa- 
rameters (f2,A, c). The equations are symmetric under the transformations {wi,qi -(->■ ^2,92}, {^jC, 91,92 ^ 
-A, -c, -91,-92}, and {n,qi,q2 -^2, -91, -92}. 

We focus on the fixed points of the system, i.e., when wi = W2 = qi ^ q2 = 0. The fixed points are given by the 
simultaneous roots of two cubic polynomials: 

fi{wi,W2) = K + l)|i + [A-c(w2 + l)]'| + ^u'i , (A5) 

f2{wi,W2) - {W2 + I)l- + [A~C{WI + 1)Y\ + —W2. (A6) 

Once wi and W2 are found, 91 and 92 can be calculated, 

9iK,W2) = r— TTT -( —7^2 ' 

I ^ -%W2[A-C{WI+1)]+I^W2 

92K,W2) = 1 -, — ■ (A8) 

Since fi(wi < —1,W2) is negative and fi{wi > 0,W2) is positive and similarly for /2 and W2, we know that the fixed 
points lie in the range wi,W2 G [—1,0]. 

By combining /i and /2, we find that the fixed points correspond to the real roots of a fifth-order polynomial h{w), 
which is too complicated to show here. Fortunately, one can factor it. Note that there are two kinds of fixed points: 
a uniform fixed point (wi — W2 and qi =92) means that the two sublattices are identical, and a nonuniform fixed 
point (wi 1x12 and 91 7^ 92) means that the system is in the antiferromagnetic phase. The uniform fixed points are 
given by real roots of a cubic polynomial, f{w) = fi{w,w), 

"rj2 1 

(A-3c)(A-c) 



f{w) = c^w^ ~ c{2A^3c)w''' 
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w + (A - c)2 + 1 . (A9) 



Since the roots of / are also roots of h, we know that / is a factor of h. Thus, h{w) = 4f{w)g{w) and the real roots 
of the quadratic polynomial g(w) correspond to nonuniform fixed points, 

g(w) = c^{l + 4:A^ + 2n^)w^-2c[{A-c){l + 4:A^) + {2A~c)n^]w 

+c^{l + 4A2) - 2cA(l + 4A2 + 2rj2) -I- i(l + AA^ + 2n^f . (AlO) 

Hence, there are at most three uniform fixed points and two nonuniform fixed points. One should think of the two 
nonuniform fixed points as being a joint pair, since they correspond to wi and W2- 

At this point, the uniform and nonuniform fixed points can be found by numerically solving for the roots of / and 
g, and their stabilities can be determined by calculating the eigenvalues of the Jacobian for each fixed point. However, 
to obtain general results, we derive as much information as possible analytically without explicitly solving for the 
fixed points. In particular, we care about the number of each kind of fixed point and their stabilities as a function of 
the parameters. 
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1. Number of uniform fixed points 

Here, we examine the number of uniform fixed points, i.e., the number of real roots of f{w). Since / is cubic, it has 
one or three real roots (two in special cases). Suppose c and A have opposite signs. The polynomial /(w) = /(w = 
w — 1) has coefficients with signs + + H — . By Descartes' rule of signs [H, f{w) has exactly one positive root, which 
means that /(w) has exactly one root with w > —1. Since we know that all roots of f{w) are in the interval [—1, 0], 
this shows that if c and A have opposite signs, there is only one root. 

Now we check when / has three roots. Since / cubic, it has three roots when the local maximum and minimum 
exist and are positive and negative, respectively. Thus, there are three roots if and only if 

4A2 > 60^+3 and (All) 



/ -2A(18172 + 4A2 + 9) - (4A2 - 6n^ - 3)i -2A(181l2 + +g) + {4A^ _ en^ _ 3)1 

According to this condition, for large |A|, there is exactly one root, i.e., one uniform fixed point, regardless of the 
sign of c. 



2. Stability of uniform fixed points 

We check the linear stability of the uniform fixed points to perturbations. Since Eqs. (jAip - (jA4l) are symmetric 
between 1 and 2, the eigenvectors of the Jacobian for uniform fixed points are either symmetric or antisymmetric 
between 1 and 2. (The symmetric eigenvectors correspond to perturbations that affect 1 and 2 identically, while the 
antisymmetric eigenvectors represent perturbations that affect 1 and 2 in opposite directions.) This is convenient, 
because we can check the stability to symmetric and antisymmetric perturbations separately, and the characteristic 
polynomials are cubic instead of sixth-degree. A uniform fixed point is stable overall if it is stable to both kinds of 
perturbations. 

The Routhe-Hurwitz criterion is very useful here because it can provide stability information without explicitly 
knowing the fixed point ^3J. Suppose the characteristic polynomial for a fixed point is cubic: a^X^ + a2X^ + aiX + ag. 
All the eigenvalues have negative real part if and only if a2,ao,aia2 — ag > 0; this means the fixed point is stable. 
All the eigenvalues have negative real part except for a pair of purely imaginary roots if and only if 02 , ag > and 
a2ai — qq = 0; this indicates a Hopf bifurcation 



a. Stability to symmetric perturbations 

First we consider the stability of a uniform fixed point to symmetric perturbations. In this case, we can consider a 
simpler system by letting w = wi — W2, q = qi = 92, and 

w = ~2nimq-w-l, (A13) 

q = i[A-ciw + l)]q-^+i^w , (A14) 

whose fixed points are given by the roots of /(w) in Eq. (jA9p . Once w is found, q is given by 



_ -fw[A-c(w + l)] 
= i + [A - c(z« + 1)]2 

The characteristic polynomial for a fixed point w,q is 



(A15) 



f 5 1 1 

a{X) = X^ + 2X^ + l[A~ c{w + 1)]^ + + - ~ 2ncRe qV X + [A - c{w + 1)]^ + — + --2ncReq. (A16) 

For a{X), we see that a2 > and a2ai — ag = oq + fl^ + 2. According to the Routhe-Hurwitz criterion given above, 
the fixed point is stable to symmetric perturbations if and only if ag > 0. Also, since a2ai — ag > ag, there is never a 
Hopf bifurcation from symmetric perturbations for any uniform fixed point. 
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Now suppose c and A have opposite signs. We showed in Sec. lA II that in this case, there is one uniform fixed 
point. We also see that 2ilcKe g < in this case and hence ag > 0. Thus the single uniform fixed point is stable to 
symmetric perturbations. 

Since / is cubic, fixed points of Eqs. (|Al3[) - (jA14|) appear and disappear through saddle- node bifurcations as the 
parameters change [ij. In a saddle- node bifurcation, two fixed points with opposite stabilities appear or disappear 
together. Also, since there is never a Hopf bifurcation, a given fixed point has the same stability as the parameters 
change. 

These statements allow us to deduce the stabilities of all the uniform fixed points. Suppose we start with c and 
A having opposite signs. There is a single uniform fixed point and it is stable as shown above. As we change the 
parameters, eventually it collides with another uniform fixed point via a saddle-node bifurcation, so the second fixed 
point must be unstable. Furthermore, that unstable fixed point undergoes a saddle-node bifurcation with a third fixed 
point, so the third fixed point must be stable, even when it is the only fixed point. (Remember that these stabilities 
are with respect to symmetric perturbations.) 

To summarize, when there is one uniform fixed point, it is stable to symmetric perturbations. When there are three 
uniform fixed points, the outer two (highest and lowest values of w) are stable to symmetric perturbations, while 
the inner one (middle value of w) is unstable to symmetric perturbations. Of course, the outer fixed points could be 
unstable to antisymmetric perturbations, which is what we check in the next section. 



b. Stability to antisymmetric perturbations 

We now check when the uniform fixed points become unstable to antisymmetric perturbations. Let the fixed point 
be w,q. We consider antisymmetric perturbations around it: wi = w + Sw, W2 — w — Sw, qi = q + Sq, and q2 = q — Sq. 
By plugging into Eqs. (jAip - (|A4p and linearizing for small Sw,Sq, we find the characteristic polynomial, 

/3(A) = + 2X^ + \ [A- c{w + 1)]2 + + - + 2ncRc q} \ + [A - c{w + 1)]^ + _- + - + 2f7c Re g , (A17) 



which is similar to Eq. (IA16|) . For /3(A), we see that 02 > and 0201 — oq = uq + H.'^ + 2. According to the 
Routhe-Hurwitz criterion given above, the fixed point is stable to antisymmetric perturbations if and only if ao > 0, 
i.e., 

02 1 

[A-c(u; + l)]2 + — + - + 2f7cReg>0. (AI8) 

We can simplify this using the fact that the fixed point satisfies f(w) — 0: the uniform fixed point is stable to 
antisymmetric perturbations if and only if 

(t){w) = c^w^ + 2c^w + c^-A^- — - -<0. (A19) 
So when a uniform fixed point is on the verge of instability, it satisfies 



2\c\ 

For large |A|, (j> < 0, so the one uniform fixed point that exists is stable to both symmetric and antisymmetric 
perturbations. 

Note that since 0201 — oq > Oq, there is never a Hopf bifurcation from antisymmetric perturbations. Since we 
already ruled out Hopf bifurcations from symmetric perturbations, we conclude that uniform fixed points never have 
Hopf bifurcations. 



3. Nonuniform fixed points 

The nonuniform fixed points are given by the roots of the quadratic polynomial giw) in Eq. (IA10[) . The two 
nonuniform fixed points appear and disappear together as the parameters change. The symmetry in Eqs. (jAl[) - (|A4[) 
between 1 and 2 means that the two nonuniform fixed points have the same stability. Thus, the nonuniform fixed 
points must appear via a pitchfork bifurcation from a uniform fixed point as the parameters change [ij. In other 
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words, the nonuniform fixed points intersect with a uniform fixed point, which changes stabihty at the intersection. 
(In Sec. lA 4[ we wih determine whether the bifurcation is supercritical or subcritical and which uniform fixed point is 
involved.) 

The intersection of the nonuniform fixed points with a uniform fixed point can be shown explicitly. From gi^w), the 
nonuniform fixed points exist if and only if A, c) < 0, where 

C(l^,A,c) = 16{l + 2n^)A'^ -32cn^A^ + 8{l + 2n^)^A'^ -8cn^{l + 2n^)A + {l + 2nY -Ac^n-^ . (A21) 

On the verge of the appearance of the nonuniform fixed points, C = ^-nd the root of g{w) is 

A 

w = -IH \ -5 ^. (A22) 

c 1 + 4A2 + 2f72 ^ ' 

One can show that Eqs. (jA20[) and (jA22[) are equal using f{w) = 0. Thus, the change in stability of a uniform fixed 
point coincides with the appearance of the nonuniform fixed points. 

In the case when |A| is large, ( ^ A**, so nonuniform fixed points do not exist. 

There is a convenient sufficient condition for the existence of the nonuniform fixed points. If ^(A = 0) = (1 + 
2172^3 _ /^(,2^4 ^ there is a range of A around A = 0, for which the nonuniform fixed points exist. In other words, 
if 



(2172 + 

there is a range of A around A = 0, for which nonuniform fixed points exist. 



c| > , (A23) 



4. Connection between uniform and nonuniform fixed points 



We describe what happens as A changes. Without loss of generality (due to symmetry), assume c > 0. Suppose A 
starts out large and negative. There is one uniform fixed point and it is stable. Nonuniform fixed points do not exist. 
Then let A increase. At some point, a uniform fixed point may undergo a pitchfork bifurcation: it becomes unstable 
to antisymmetric perturbations and the nonuniform fixed points appear. The fact that the nonuniform fixed points 
exist when the uniform fixed point is unstable indicates that the bifurcation is supercritical [l|. 

When the pitchfork bifurcation happens, there may be three uniform fixed points. Which one undergoes the 
bifurcation? Since (p{w) ~ in Eq. (jAlQp . as A decreases, the first uniform fixed point to go unstable must be 
the upper one (the one with the highest w). Furthermore, since there can be at most two nonuniform fixed points, 
only the upper uniform fixed point may be unstable to antisymmetric perturbations. We conclude that when there 
are three uniform fixed points, only the upper one may undergo a pitchfork bifurcation, and the lower one is stable. 
According to Sec. lA 2 al the middle uniform fixed point is always unstable. 

It is possible that the pitchfork bifurcation happens when there is only one uniform fixed point. Then obviously 
that fixed point must undergo the pitchfork bifurcation. 

Since the pitchfork bifurcation is supercritical, the nonuniform fixed points are stable when they appear. 

Then let A increase towards 00. Eventually C > and (f> < again. This indicates that there was another pitchfork 
bifurcation in which the nonuniform fixed points disappeared and the uniform fixed point became stable again. 

Note that we have not determined here whether the nonuniform fixed points ever have Hopf bifurcations. As stated 
in the main text, we found numerically that they do sometimes have Hopf bifurcations. 



Appendix B: Complete lattice 



Here, we study the solutions for the complete lattice, keeping track of every site instead of lumping them into 
sublattices. We show that there is a dynamical instability, in which the uniform steady state becomes unstable to 
perturbations of wavelength 2. The solutions of the complete lattice are the same as the solutions using sublattices, 
so the antiferromagnetic transition is not an artifact of using a bipartite lattice. 

The approach is to find the uniform steady state and see when it becomes unstable to perturbations. Consider a 
d-dimensional lattice with N sites in each direction. Let n be a d-dimensional position vector. The system is described 

by 



Wfi 



-2f2 Im (7,-5 



1 



A-bJ2 + 1) 



(rhn) 



2 



n 



(Bl) 
(B2) 
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where h = V/2^ = c/2d\s the nearest- neighbor interaction. 

The uniform steady state {wfi = w, qfi = q) is given by the fixed points of the system, 

w = -2f7Im q-w-l , (B3) 

q n 

q = i [A - 2db{w + 1)] g - I + i^w . (B4) 

These equations are the same as Eqs. (|Al3P - (jA14|) but with c = 2dh. Thus, we can use the resuhs of Sec. lA 2 al The 
uniform steady state is given by the real roots of /(w) in Eq. (|A9|) . and there are one or three solutions. The Jacobian 
of Eqs. (jB3P - (|B4[) determines the stability of a uniform solution to uniform perturbations, i.e., an identical offset to 
every site. 

Now we consider perturbations 5wfi, 5qfi around the uniform steady state: 

Wfi ^ w + dwfi , (B5) 
9n = q + Sqn ■ (B6) 

We write them in terms of Fourier components Sw^, Sqj:: 

k 
k 

k = (fci,fc2,...,fcd) , (B9) 
ki = ^J, j-0,...,Ar_l . (BIO) 

We write Eqs. (|B1I) - (|B2I) in terms of S'Wj:,Sqj: and linearize. The Fourier components are uncoupled from each other 
and each component k is described by three differential equations: 

S^j: = - Srj:) - Swj: , (Bll) 

% = -2i5g^^cos/c<?^ (5i&^r + i^5wj + |i[A-2d6(w + l)]-i|%, (B12) 

dq*_^ = 2i6g* ^^cosfcf^ + |-i[A- 2d6(w + 1)] - ij^g*- . (B13) 

The characteristic polynomial for component k is 

ri{X) = + 2X'^ + i^[A- 2db{w+ l)f + + ^-AnbRe qJ2 cos ki^ X (B14) 

+ [A-2db{w + l)]^ + — + --4nbRe q'^coski . (B15) 



Note the similarity to Eqs. (IA16P and (|A17I) . Now we use the Routhe-Hurwitz criterion given in Sec. lA 21 For r]{X), 
a2 > and 0201 — Uq = Qq + il^ + 2. Thus, the uniform steady state is stable to a perturbation with wave vector k if 
and only if 



[A - 2db{w + 1)]^ + — + - - AnbRe q'^coskg > . (B16) 



Suppose the uniform steady state satisfies 4i76Re q > 0. As the parameters change, the first mode to go unstable is 
the one with all ki — 0. This corresponds to uniform perturbations that simply offset the uniform solution. Thus we 
identify this uniform steady state as the unstable fixed point of Eqs. (IB3P - (jB4p . Based on the discussion in Sec. lA 2 al 
this uniform solution is actually always unstable when it exists. 
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Now suppose the uniform steady state satisfies 4f26Re q < 0. As the parameters change, the first mode to go 
unstable is the one with all ki = tt. This corresponds to perturbations of wavelength 2, i.e., antifcrromagnetic 
perturbations. The uniform steady state is unstable to this mode when 



This is the same as Eq. (|A18P with c = 2db. Thus, the discussion in Sec. lA 41 applies here: as the parameters change, 
the uniform solution becomes unstable to the antiferromagnetic solution. One can find the antifcrromagnetic solution, 
but that is equivalent to the nonuniform fixed points of Eqs. (jAl[) - (jA4[) . i.e., the real roots of g{w) in Eq. (IA10[) . 
Thus, the complete lattice has the same antiferromagnetic transition as the bipartite lattice. 
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